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ABSTRACT 

The broad Iron line in MCG-6-30-15 is fitted to the Comptonization model 
where line broadening occurs due to Compton down-scattering in a highly 
ionized optically thick cloud. These results are compared to the disk hne model 
where the broadening is due to Gravitational/Doppler effects in the vicinity of 
a black hole. We find that both models fit the data well and it is not possible 
to differentiate between them by fitting only the ASCA data. The best fit 
temperature and optical depth of the cloud are found to be kT — 0.54 keV and 
r = 4.0 from the Comptonization model. This model further suggests that while 
the temperature can be assumed to be constant, the optical depth varies during 
the observation period. We emphasis an earlier conclusion that simultaneous 
broad band data (3 — 50 keV) can rule out (or confirm) the Comptonization 
model. 

Subject headings: accretion disks — black hole physics — galaxies individual 
(MCG-6-30-15) — galaxies : Sey f ert — line :profile 
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1. Introduction 

A long (rs 4.5 days) observation by ASCA of the Seyfert 1 MCG-6-30-15 revealed that 
the Iron hne profile in this source is broad with velocity width of > 10^° cm sec~^(Tanaka 
et al. 1995). A more detailed analysis of the data by Iwasawa et al. (1996) showed that the 
profile is variable and the line width is maximum when the source intensity is minimum. 
Broad Iron lines have also been detected in other AGN by ASCA ( e.g. Nandra et al. 1997). 
The most straight forward interpretation of this feature is the disk line model where the 
line is broadened due to the combined effects of Doppler and gravitational red-shift in the 
vicinity of a black hole (Fabian et al. 1989 ). Prom modeling the X-ray data, Iwasawa et al. 

(1996) showed that the disk producing this line has to have an inner radius 1.2rs where 
Tg = GM/c^ is the Schwarzschild radius. This would imply that the black hole is spinning 
close to its maximal value. If this interpretation is correct then this would be the first 
direct observation of the strong gravitational effects expected in the vicinity of a black hole. 
This would strongly constrain theoretical models since the inner disk region has to be cold 
in order to produce the Iron line while at the same time the hard X-ray producing region 
would also have to be located close to the inner edge of the disk. Reynolds & Begelman 

(1997) have argued that if some of the Iron hne emission arises from the region inside the 
last stable orbit, the black hole need not be a rotating one. 

Sulentic et al. (1998) showed there was a disagreement between the inclination angle 
derived by fitting the Iron line profile and that from HI and Ha measurements. Several 
AGN for which a broad Iron line was detected have the peak centroid close to 6.4 keV 
(Nandra et al. 1997). It was pointed out by Sulentic, Marziani & Calvani (1998) that 
this is not expected if the disks are oriented in random directions. They propose that 
the Iron line profile is a sum of two independent Gaussian lines in order to explain the 
different correlation relationships of the red and blue side with the source intensity. They 
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also claimed that there is a blue wing in the Iron line emission which cannot be explained 
by the disk line model. These discrepancies and the importance of the implication of the 
disk line model warrant a study of alternate models to explain the phenomenon. Due to the 
profound implications of being able to 'observe' the immediate vicinity of a supermassive 
black hole, it is important to examine models for the broad Iron line that do not require 
strong gravity. 

As an alternative to the disk hne model, Czerny, Zbyszewska and Raine (1991) 
proposed that the line is intrinsically narrow and it gets broadened due to Compton down 
scattering of the photons as they pass through an optically thick cloud. This model is 
referred to here as the Comptonization model. Both, the Comptonization and the disk line 
models, predict a broad Iron line with an extended red wing. However, the detailed spectral 
shape of the line in the two models is different. In particular the disk line model predicts a 
double peaked line profile while the Comptonization model predicts a smooth feature. The 
observed line profile in MCG-6-30-15 can be fitted well by the disk hne model (Iwasawa et 
al. 1996), but so far, no fits of this data have been made to the Comptonization model. 
Fitting the data to the Comptonization model would yield values for the optical depth and 
the temperature of the cloud - parameters that can be used to test the physical feasibility 
for the existence of such a cloud. 

Fabian et al. (1995) rejected the Comptonization model by arguing that the 
surrounding Comptonizing cloud has to have a radius R < 10^^ cms in order that the cloud 
be highly ionized and does not produce strong absorption lines. For a IO^Mq black hole 
this would imply that the Iron line producing region is smaller than 50rs and gravitational 
effects would be important. The lack of a blue wing in the best disk-line fit to the hne 
profile, implies that the temperature of the cloud should be kT < 0.2 keV, which is in 
apparent confiict with the expected Compton temperature of the cloud. Further a break 
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in the continuum around 20 keV is expected for the Comptonization model which has 
not been observed. On the other hand, Misra & Kembhavi (1998) pointed out that for a 
smaller sized black hole, the intrinsic Iron line produced may not be significantly broadened 
by gravitational effects. They calculated the equihbrium temperature of the cloud to be 
around 0.2 keV provided an intense UV source is assumed to be present in this source. 
They showed that the present broad band data for this source is consistent with the 
Comptonization model. In these arguments the temperature of the cloud was constrained 
by the absence of a blue wing in the fitted spectrum for the disk line model. However, the 
unfolded spectrum depends on the model used to fit the data. In other words, the inferred 
shape of the Iron line profile depends on the model used to describe the fit. It is important 
to obtain better estimates for the temperature and optical depths of the cloud in order to 
verify the arguments presented by Fabian et al. (1995). In this paper, we directly fit the 
observed line profile to the Comptonization model and hence obtain more realistic values of 
the cloud parameters. 

The Comptonizing model raises several questions regarding the stability, geometry and 
dynamics of such a cloud. The origin of such cloud(s) is unspecified but some arguments 
for their existence are presented by Guilbert & Rees (1988) and Kuncic, Blackman & 
Rees (1996). Optically thick obscuring clouds may be responsible for the intensity dip in 
MCG-6-30-15 (McKernan & Yaqoob 1998 ; Weaver & Yaqoob 1998). 

The broad line may also be an artifact of simplistic modeling of the continuum. If the 
underlying continuum is complex and not a simple power-law, residuals caused by this fit 
may have been interpreted as a broad line. In order to test this hypothesis, we fit the data 
to a narrow line and a broken power-law as the continuum. 

In the next section we briefiy describe the observations and data analysis technique. In 
§3 the three models used for the fitting are described. In §4 the results of the data analysis 
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is presented. In §5 we discuss the implications of the results and summarize the paper in §6. 

2. Observations 

ASCA observed MCG-6-30-15 from 1994 July 23 to 27 (Tanaka, Inoue & Holt 1994) 
and the data has been analyzed in detail by Iwasawa ct al. (1996). The SIS data which has 
been analysed in this paper, was obtained in Faint/ICCD mode. We used screened events 
files supplied by the HEASARC onhne service for this observation, where the Faint mode 
data was already converted into Bright (B) and Bright2 (B2) modes for analysis. In this 
paper, we analyze the 3 — 10 keV data since below 3 keV the spectrum is affected by the 
partially ionized gas ("the warm absorber") surrounding the source. Following Iwasawa et 
al. (1996), the data set was divided into three epochs of which the epoch labeled "high" 
encompasses the highest intensity level of the source. The epoch labeled "low" encompasses 
the lowest intensity level of the source. The rest of the observation was grouped into a single 
data set labeled "medium" . Data from both the SIS chips (SIS 1 and 2) for the Bright 
and Bright2 modes were grouped and analyzed together. Only the relative normalization 
between these four sets of data was allowed to vary, but in all cases the variation was found 
to be less than 2%. 

3. Brief description of Models 

3.1. Disk Line model 

The disk line model of Laor (1991) assumes that the emission arises from an accretion 
disk around a black hole and includes general relativistic effects. The inner radius of the 
accretion disk is rj, where = 1.235rs and = 6rs (r^ — GM/c^) correspond to the 
case of a maximally rotating black hole and a non-rotating one respectively. The model is 
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parameterized by r^, the outer radius Tq, the emissivity index ^, the rest frame hne energy 
(fixed at 6.4 keV) and hne intensity. The emissivity index ({) specifies the dependence of 
the emissivity with radius. The continuum is assumed to be a power-law characterized by 
two parameters - photon index F and a normahzation factor. Thus the total number of free 
parameters used for the fitting is six. 

3.2. Comptonization model 

The physical geometry of this model is that the central engine consists of a hard 
X-ray source and a cold (UV producing) medium. The X-rays impinge upon the cold 
medium and produce a fiuorescence Iron line and a refiected spectrum. However, unlike 
the disk line model, the Iron line is produced sufficiently away from the black hole and is 
intrinsically narrow and is broadened by the presence of an extended Comptonizing cloud. 
It is important to note that for self consistency both the intrinsic power-law and Iron line 
are Comptonized by the surrounding cloud. In principle, a reflected component which 
peaks around 50 keV should also be taken into account. Since the ASCA observations are 
restricted to photons with energy less than 10 keV we can neglect the reflected component 
here (Iwasawa et al. 1996). 

The intrinsic spectrum is assumed to be a power-law with photon index F and a narrow 
Gaussian hne with a — 0.05 keV and centroid energy E — 6.4 keV. The source is assumed to 
be surrounded by a highly ionized cloud with optical depth r = narR at a temperature kT. 
The output spectrum is calculated using the Kompane'ets equation. The Comptonization 
model has seven parameters - r, kT, the intrinsic photon index F, the line energy (flxed 
at 6.4 keV ), the intrinsic width of the hne (flxed at 0.05 keV), the line intensity and a 
normalization factor. Thus the total number of free parameters used for the fltting is flve. 



-8- 



3.3. Broken power-law model 

The broken power-law is an empirical model used to test the hypothesis that the Iron 
line in this source is narrow and the artificial broadening observed is due to complexities in 
the underlying continuum. The Iron line is assumed to be a Gaussian with width a = 0.2 
keV and centroid at 6.4 keV. The continuum is fit by a broken power- law characterized 
by a high energy photon index Fi, a low energy photon index a break energy Eb and 
normalization. Thus the total number of free parameters used for the fitting is five. 

4. Results 

The medium intensity (MI) data were fitted to the disk line model of Laor (1991) for 
later comparison with the Comptonization model. The best fit parameters to the former 
model are shown in Table 1 and the unfolded data is plotted in Figure 1. As mentioned 
by Iwasawa ct al. (1996) this model fits the data well with xV (dof)= 1377/(1421) (Table 
1, first row). The inner radius ri{— 2.9rs) is constrained to be less than 3.7rs. However, 
since the line is broader during the low intensity phase, the inner radius is probably close 
to 1.235 which corresponds to a maximally rotating black hole. Setting = 1.235rs and 
refiting the data gave Ax^ = 1-0 (Table 1, second row). 

The low (LI) and high intensity (HI) data sets were fitted to the disk line model after 
setting Ti = 1.235, = 16.4 and inclination angle 4> = 32°. These parameters are not 
expected to change during the observation period and freezing them constrains the rest of 
the parameters. Moreover the the number of free parameters become equal to those of the 
Comptonization model (see below) which enables a direct comparison . The HI data set is 
formally well fit by the model with x^/ (dof)= 1235/(1424) (see Table 1). There seems to 
be evidence for an additional narrow component (Iwasawa et al. 1996 ; Sulentic, Marziani 
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& Calvani 1998). The emissivity index = 2.0l'lo;7) does not vary significantly from the MI 
state. Much of the spectral change can be attributed to the slight change in the power-law 
slope. The line shape during the LI level is different from that of the MI data (Table 1). 
The hne is broader as indicated by the increase in the the emissivity index = S.Sioie); 
since an higher ^ imphes that more photons are produced closer to the Black hole. The 
data is again well fit by the model with x^/ (dof)= 967/(1424). 

The Comptonization model also fits the MI data well with xV (dof)= 1396/(1423) 
although with a higher x^ than the disk line model (Table 2). The temperature of the 
Comptonizing cloud is found to be kT — 0.541q;^9 keV. This is hotter than the temperature 
estimated by best fit hne profile for the disk hne model which was kT < 0.2 keV (Fabian 
et al. 1995, Misra & Kembhavi 1998). The reason for this discrepancy is that fitting the 
data to the Comptonization model reveals a blue wing (Figure 2). The inferred spectral 
shape i.e. the unfolded spectra depends upon the model used to fit the data. Thus the 
absence of a blue wing in the unfolded spectrum for the disk line model cannot be used as 
an estimation of the temperature of the cloud. The LI and HI data sets were fitted after 
setting kT = 0.54. The HI data is well fit by this model xV (dof)= 1225/(1424) (see Table 
2). The fit is similar to the disk hne model and again there seems to be evidence for an 
additional narrow component (Figure 3). Compared to the MI fit the optical depth and line 
strength have increased while the power-law index (F) remain unchanged within the limits 
of error. Thus in terms of the Comptonization model, the photon index variation in the fit 
to the disk line model can be attributed to the change in degree of Comptonization of the 
intrinsic continuum. However, the limited statistics of the data set and the presence of an 
additional narrow line does not allow for any concrete statements to be made. The LI is 
well fit with a larger optical depth r = 6.7 with x^/ (dof)= 969/(1424). The larger optical 
depth leads to a broader line for this set. 
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Table 1: Spectral parameters for the disk line model of Laor (1991) for the intensity selected 
data sets. The parameters are the power- law photon index (F), the emissivity index (^), the 
inner radius (r,) in r^, the outer radius (ro) in r^, the line Intensity (/ane in units of 10~^ 
photons/sec/cm^ and the inclination angle in degrees. The line energy is fixed at 6.4 keV. 
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Table 2: Spectral parameters for the Comptonization model for the intensity selected 

data sets. The parameters are the power-law photon index (F), the optical depth (r), 

the temperature of the cloud {kT) in keV and the line Intensity {lune) in units of 10~^ 

photons/sec/cm^. The line energy is fixed at 6.4 keV and the intrinsic line width is assumed 

to be cr = 0.05 keV. The difference in between the disk fine (Table 1) and Comptonization 

model is tabulated in the last column. 

Data set T t kT lu^e xV (dof) Ax^ 

Medium imtf^^ 4.16lg:l^ 0.541^:^9 2.221^:^ 1396/(1423) +18 

High 1.96l:Ji 5.32lg:^^ 0.54 3.711^:^^ 1225/(1424) -10 

Low 1.761:^1 6.691?:^^ 0.54 3.8ll?:^^ 969/(1424) +2 
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Since the data is from a narrow band (3 - 10 keV) both spectral slopes of the broken 
power-law model cannot be constrained. Therefore, the MI data is fit to the broken 
power-law model with the high energy photon index fixed at 2.0. This gives an unacceptable 
fit with I (dof)= 1506/(1424). An Acceptable fit was obtained when the high energy 
photon index was set to 2.5 with x^/ (dof)= 1422/(1424). The HI and LI data sets can also 
be fitted by this model provided the high energy photon index is 2.5 (Table 3). 



5. Discussion 

The fit to the disk line model to the MI and LI data set is marginally better than 
for the Comptonization model. The reduced for both the models are less than unity. 
The reason for this degeneracy is that the spectral shape of the best fit disk line model 
does not have a clear double peaked feature (Figure 1). Instead the red wing of the line 
gradually extends to w 5 keV and such a feature can also be due Compton down-scattering 
of photons. Thus, the narrow band (3 - 10 keV) data used here cannot differentiate between 
these two models. 

Although the HI data set is formally well fit by the models, there is evidence for an 
additional narrow Iron line (Iwasawa et al. 1996) in the residual plots. This could be due 
to a real additional narrow Iron line ( Sulentic, Marziani & Calvani 1998) or because the 
Iron line shape is variable in a time scale shorter than the observation time 10^ sees). In 
the latter scenario the observed shape of the Iron line would be the time-averaged profile 
and therefore could not be described by time-independent modeling. 

Since there are no strong X-ray absorption features in this source the size of the 
Comptonizing cloud has to be smaller than R < 10^^ cms ( Fabian et al. 1995). For a 
IO^Mq black hole this corresponds to about ~ 300r5. The intrinsic Iron line produced in 
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this case will not be significantly broadened by gravitational red-shift or Doppler effects. 
However if the black hole mass is larger, these effects would be important and the observed 
profile would be complex combination of Comptonization and gravitational effects. In this 
work for simplicity we have assumed that the intrinsic line is narrow {a — 0.05 keV). The 
temperature of this cloud can be determined by equating the heating and cooling of the gas 
by Compton scattering. As described by Misra & Kembhavi (1998) the temperature of the 
gas depends upon the assumed UV flux of the intrinsic source. In Figure 4, the dashed line 
shows the assumed intrinsic spectrum for the source. The solid line is the output spectrum 
after the radiation passes through the cloud with optical depth r = 4.0 (the best fit value 
for the medium intensity data) . The temperature of the cloud was calculated by balancing 
the input radiative power to the output power. The UV flux has been assumed here such 
that this temperature is equal to the best fit value kT 0.55 keV. The UV luminosity 
is lower than that estimated by Misra & Kembhavi (1998) since the best fit value of the 
cloud temperature for the Comptonization model is 0.5 keV and not 0.2 keV as has been 
assumed earlier. The shape of the UV bump has been assumed here to have a black-body 
like shape. However, the equilibrium temperature is not sensitive to the spectral shape but 
rather to the total luminosity in the UV band compared to the X-ray band. The UV flux 
in MCG-6-30-15 is highly reddened and hence a UV bump is not directly observed for this 
source (Reynolds et al. 1997). Using the minimum value for the extinction, a lower limit 
of the intrinsic flux (F^^n) has been calculated by Reynolds et al (1997) and is plotted in 
Figure 4 for comparison. The UV flux required by this model is about a factor ten larger 
than this minimum value. The reprocessed IR emission from covering dust in this source 
also indicate that the UV flux is much larger than Fmin (Reynolds et al 1997). 

The temperature of the gas will react to the changing intrinsic spectrum in a time scale 
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set by the time for the photons to diffuse out of the system, 



tc ^ (— = 8 X 10Ws( 




(1) 



The dynamical time scale is 
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where v « {2GM/R) 



is the average bulk velocity of the plasma. Since this two time-scales 



are comparable to the observation time scale , the temperature and the optical depth can 
vary during the observations. The data is compatible with no change in temperature during 
the observation period with the change in the line profile being attributed to changes in the 
optical depth of the cloud ( Table 2). 

The continuum spectrum for energies {E > 20 keV) is greatly affected by the presence 
of a Comptonizing medium. For £■ > 100 keV the decrease in flux is nearly half of what 
is expected from a simple power-law model ( Figure 4). The averaged Seyfert 1 spectrum 
from OSSE has been presented by Gondek et al. (1996). They find that at ii^ = 100 
keV, there is no significant decrease in the flux from the extrapolated power-law value. 
However, the statistics is not good enough to rule out a 50% decrease in flux. Moreover, the 
averaged spectrum may not represent the actual spectrum of MCG-6-30-15 at a given time. 
Significant differences in the continuum spectrum is also expected for 20 < £^ < 100 keV. 
The analysis in this energy range is complicated by the presence of the reflection component 
which is significant only for 15 < £^ < 60 keV. Nandra & Pounds (1994) analyzed the 
Ginga data for this source in 2 < < 18 keV range and find that a reflection component 
corresponding to a solid angle Qj. ~ 2.07r is required to explain the data. Recently, Lee 
et al. (1998) analyzed the RXTE (PCA) data for this source for the energy band (2-20 
keV). They find a refiection component corresponding to a larger solid angle Vt^ ~ 3.07r. 
Preliminary analysis of the BeppoSax observation for this source ( Matt 1998) requires a 
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Qr ~ 2.67r. These results indicate that there is an excess flux at £■ > 15 keV which is in 
apparent contradiction to the Comptonization model. This excess can be explained within 
the framework of the Comptonization model if the intrinsic reflection component has a 
higher value of ^ 1.5 times the fitted Qr (Misra & Kembhavi 1998). For E ^ 50 keV, 
the difference between the Comptonized power-law and power-law model is larger and can 
only be removed by invoking an extremely large unphysical solid angle for the reflected 
component fir > 57r. MCG-6-30-15 has been observed by the BeppoSax PDC instrument 
which is sensitive in the energy range (13 < E < 200 keV ). Matt (1998) show the count 
spectrum for this source. However, the count flux at £^ fa 50 keV is uncertain by about 
a factor of two. Hence it may not be possible to rule out the Comptonization model by 
this observation by BeppoSax. Similar longer duration broad-band (3 < < 100 keV ) 
coverage of this source is required to rule out or confirm the Comptonization model. Such 
an observation is also possible by HEXTE instrument in RXTE. An additional complication 
arises because the source (and hence the Comptonizing cloud) is known to be variable. 
Hence simultaneous observation in this range will be required. 

Another explanation for these observations could be that the intrinsic line is narrow 
but the underlying continuum is complex i.e. not just a simple power-law. To test this 
hypothesis we have fit the data to a narrow line and a broken power-law. We find that 
the high energy spectral index has to be ~ 2.5 for this model to fit the data. The recent 
BeppoSax data (Matt 1998) seems to contradict this result since the spectral index obtained 
was 2.0. Moreover, it is also not clear which physical process will give a break energy 
^ 5 keV (Table 3). 
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6. Summciry 

The disk line model and the Comptonization model can fit the Iron profile of this 
source for the low and medium intensity data sets. The width and skewness of the Iron line 
profile can be reproduced if the source is surrounded by a cloud with optical thickness r « 4 
and temperature kT ^ 0.6 keV. This temperature is higher than earher estimates based on 
fits using the disk line model. The cloud will have this temperature provided there is a UV 
bump in the source which is not visible due to extinction. 

The Iron line profile could also be narrow if the underlying continuum is complex. 
Modeling the continuum as a broken power-law allows for a reasonable fit only if the 
high energy photon spectral is ~ 2.5 which may be in conflict with the recent BeppoSax 
observations. 

As pointed out earlier by Misra & Kembhavi (1998), we emphasis the need for 
simultaneous long duration broad band (3 — 50 keV) observations of this source. Such 
observations will be able to rule out (or confirm) the Comptonization model and will 
strongly constrain any theoretical accretion disk models invoked to explain the AGN 
phenomena. 

The authors would like to thank S. Raychoudhary and A. Kembhavi for useful 
discussions. This research has made use of data obtained from the High Energy 
Astrophysics Science Archive Research Center (HEASARC), provided by NASA's Goddard 
Space Fhght Center. 
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Table 3: Spectral parameters for the Broken power- law model for the intensity selected data 

sets. The parameters are the high energy power-law photon index (Fi), the low energy 

power-law photon index (r2), the break energy (Eb) in keV and the line Intensity {Inne in 

units of 10~^ photons/sec/cm^. The line energy is fixed at 6.4 keV and the intrinsic line 

width is assumed to be cr = 0.2 keV. The difference in between the disk line (Table 1) 

and the Broken power-law model is tabulated in the last column. 
Data set lu^, xV (dof) Ax^ 

Medium 2.5 l.Trt^ij 5.8i;[|i^ OASto^l 1422/(1424) +44 

High 2.5 IMtom 5.281^;?^ O-^Q-oH 1208/(1424) -27 

Low 2.5 lABtoit 5.30tg:i 0.27lg:i? 972/(1424) +5 
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Fig. 1. — The unfolded spectra and best fit disk line model of Laor (1991). The medium 
intensity data from both chips (SISO and SISl) for Bright (B) and Bright2 (B2) modes are 
shown separately. The dotted line is the power-law while the solid line is the sum of the 
power-law and Iron line profile. The parameters of the fit are given in Table 1. 

Fig. 2. — The unfolded spectra and best fit Comptonization model. The medium intensity 
data from both chips for Bright (B) and Bright2 (B2) modes are shown separately. The 
dotted line is the Comptonized power-law while the solid line is the sum of the power-law 
and Iron line profile. The parameters of the fit are given in Table 2. 

Fig. 3. — The unfolded spectra and best fit Comptonization model. The High intensity 
(I3)data from both chips for Bright (B) and Bright2 (B2) modes are shown separately. The 
dotted line is the Comptonized power-law while the solid line is the sum of the power-law 
and Iron line profile. The parameters of the fit are given in Table 2. 

Fig. 4. — The Comptonized multi-wavelength spectrum for MCG-6-30-15. The dotted 
line is the assumed intrinsic spectrum of the source while the solid line is the Comptonized 
one. The intrinsic component consists of a UV bump, a X-ray power-law and a refiection 
component. The fiux of the UV bump is chosen such that the equilibrium temperature of 
the surrounding cloud is kT = 0.55 keV. Also shown is the estimated UV lower hmit for this 
source ( Reynolds et al. 1997). Maximum deviation occurs for E > keV. 
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